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ABSTRACT 


Adsorbed  atoms  on  the  surface  of  a  harmonic  lattice  are  immersed  in  a  strong 
laser  field.  The  optical  Bloch  equations  are  derived,  which  include  the 
thermal  relaxation  and  the  coherent  excitation  of  the  adbond.  This  is 
accomplished  by  a  transformation  to  dressed  states,  which  diagonalizes  the 
interaction  with  the  laser.  The  single-phonon  couplings  are  then  understood 
as  transitions  between  dressed  states.  The  radiative  contributions  for 
arbitrary  strong  fields  are  obtained  in  the  master  equation,  and  it  is  shown 
that  the  coherences  with  respect  to  the  dressed  states  decay  exponentially, 
due  to  the  phonon  relaxation.  General  properties  of  the  competing  phonon- 
induced  redistribution  and  optical  excitation  of  the  level  populations  are 
presented,  and  exemplified  by  an  explicit  elaboration  of  a  three-level 
system.  The  results  are  amenable  to  analytical  evaluation  once  the 
interaction  potential  is  prescribed,  and  extensions  of  the  approach  to 
include  multiphonon  processes  are  straightforward. 
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I .  INTRODUCTION 


Relaxation  or  desorption  of  an  atom  in  a  vibrational  bond  with  a 
surface  of  a  crystal  with  finite  temperature  is  fairly  described  by  a  Pauli- 
type  master  equation  and  has  been  studied  extensively.^  ^  The  vibrational 
levels  of  the  adbond  are  coupled  by  the  phonon  field  of  the  crystal,  and 
phonon-exchange  reactions  of  the  bond  with  the  crystal,  which  acts  as  a 
thermal  bath,  cause  the  relaxation  or  desorption.  The  latter  process  can  be 
either  a  result  of  successive  resonant  single-phonon  excitations  up  the 
ladder  of  vibrational  states  or,  when  the  level  separations  are  larger  than 
the  Debye  frequency  as  a  one-step  multiphonon  process.  Besides  the 
academic  interest  in  these  processes,  the  thermal  desorption  of  adatoms  is 
the  major  technical  method  to  obtain  clean  surfaces.  However,  when  the 
distances  between  the  lower  levels  are  large  in  comparison  with  tu^,  the 
desorption  rate  might  be  very  low,  since  the  first  excitations  can  only 
occur  through  multiphonon  processes. 

An  obvious  way  to  enhance  the  desorption  is  by  irradiation  of  the 
surface  with  infrared  laser  light.  The  laser,  with  frequency  u^,  is  tuned 
in  resonance  with  a  vibrational  transition.  In  this  fashion  one  can 
efficiently  populate  a  high-lying  state,  such  that  the  transition  to  the 
continuum  can  be  established  by  a  resonant  one-phonon  process. ^  ^  There  is 
obviously  a  competing  effect,  which  diminishes  the  efficiency  of  this 
procedure.  An  optically-excited  bond  can  decay  to  a  lower  state  under 
emission  of  a  phonon  into  the  crystal.  This  can  be  considered  as  a 
conversion  of  a  photon  into  a  phonon,  and  the  net  result  is  an  energy  flux 
into  the  crystal,  without  desorption  of  the  atom.  On  a  much  larger  time- 
scale  however,  this  process  heats  the  crystal,  and  consequently  the 
desorption  is  again  enhanced. 


The  problem  of  irradiation  of  an  adsorbed  atom,  in  combination  with 

thermal  relaxation  due  to  the  coupling  with  the  phonon  reservoir,  has  been 
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treated  in  the  weak-field  limit  with  a  perturbative  approach  and  for 
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special  choices  of  the  coupling  potential  (harmonic,  Morse).  A  systematic 
development  of  a  concise  theory  for  strong  laser  fields,  which  does  not  rely 
on  a  specific  form  of  the  potential,  is  apparently  not  available  yet.  In 
this  paper,  we  investigate  the  interaction  of  an  intense  radiation  field 
with  the  bounded  atom.  We  interpret  the  results  in  terms  of  transition 
diagrams,  and  we  identify  the  contributions  of  pure-phonon,  phonon-photon 
exchange  and  one-phonon/two-photon  processes  to  the  laser-assisted 
redistribution  of  the  level  populations  and  to  the  decay  of  the  coherences. 
Common  features  are  illustrated  by  an  example. 

II.  EQUATION  OF  MOTION 

We  consider  an  atom  with  mass  m,  which  is  adsorbed  on  the  surface  of  a 
harmonic-lattice  crystal.  The  van  der  Waals  forces  are  represented  by  the 
potential  V  between  the  adatom  and  its  nearest  surface  atom,  and  by  an 
effective  surface-potential  $.  The  lateral  motion  of  the  atom  will  be 
neglected,  since  it  couples  only  weakly  to  the  phonon  field,  and  the 
direction  perpendicular  to  the  surface  will  be  denoted  by  z.  If  we  choose 
the  origin  of  our  coordinate  system  at  the  equilibrium  position  of  the 
surface  atom,  with  mass  M,  and  indicate  the  position  of  m  by  zej_  and  of  M  by 
u,  then  V  will  only  depend  on  the  distance  |zej_  -  u|.  Since  |u|  is  much 
smaller  than  |z|,  we  can  expand  the  potential  as 

V( | ze^  -  u|)  =  Vq( z )  -  u.e^  —  VQ(z)  *  (2.1) 

where  Vq(z)  is  defined  as  the  interaction  for  u  =  0.  We  will  only  retain 
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riz 


these  two  terms,  which  implies  that  we  neglect  multiphonon  processes.  This 
will  keep  the  formulation  concise  and  the  interpretation  transparent.  It 
is,  however,  straightforward  to  retain  higher-order  terms.  The  atomic  part 
of  the  interaction  is  abbreviated  as 


s-eV>  • 


(2.2) 


and  the  amplitude  of  the  phonon  field  for  the  surface  atom  M  will  be 


represented  by 


Kfiv1  M/2  (  tv 

2Mvto(k) J  aks  aks  ~ 


ks 


(2.3) 


Here  v  and  v*  are  the  volumes  of  the  crystal  and  unit  cell,  respectively, 
e^s  is  the  unit  polarization  vector  of  the  phonon  mode  ks,  a^s  and  a*s  are 
the  annihilation  and  creation  operator  for  the  mode  ks,  and  u>(k)  is  the 
dispersion  relation.  With  this  notation,  the  vibrational  coupling  of  the 
adbond  to  the  crystal  is  accounted  for  by  the  Hamiltonian  K  +  +  $  +  H  ^ 

u.ejS,  with  K  the  kinetic  energy  of  m  and  H  ^  the  free-field  phonon 


Hamiltonian , 


V  'I 


'f“<k)  “ks  aks 


(2.4) 


The  pure  atomic  part  K  +  +  $  of  the  Hamiltonian  will  be  termed  the 

bare-atom  Hamiltonian,  since  it  describes  the  atomic  bond  with  u  =  0  and 
without  a  laser  field.  Its  eigenvalues  -fiuK  and  eigenstates  |i>  are  defined 


(K  +  V  +  $)|i>  =  *u>.  |  i>  , 


(2.5) 


and  explicit  representations  in  coordinate  space  can  easily  be  found  for  a 
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number  of  potentials.  A  prime  example  here  is  the  Morse  potential,  which 
provides  a  fair  representation  of  the  binding  potential  well,  and  for  which 
analytical  expressions  for  -fiuK  and  |i>  are  known.  With  respect  to  its  own 
eigenvectors,  we  can  write  the  bare-atom  Hamiltonian  as 


K  +  v„  + 


ipi  • 


(2.6) 


P.  =  |  i><i I  . 


(2.7) 


the  projector  onto  the  eigenstate  |i>.  We  shall  assume  that  for  the  atomic 
part  S  of  the  interaction  with  the  phonon  field,  the  diagonal  matrix 
elements  <i|s|i>  vanish.  This  is  only  an  approximation  in  the  presence  of  a 
laser  field,  since  the  <i|S|i>-term  will  give  rise  to  a  small  radiative 
contribution,  but  not  to  pure  phonon  transitions. 

A  continuous  single-mode  laser  is  tuned  into  resonance  with  a  specific 
transition  of  the  vibrational  bond.  The  electric  component  of  the 
electromagnetic  field  at  the  position  of  the  bond  is 


E(t)  =  Eq  Re[e  exp(  -  ioo^t )  ] 


(2.8) 


and  the  coupled  levels  will  be  denoted  by  |e>  and  |g>.  This  excited  state 

and  ground  state  (which  is  not  necessarily  the  lowest  state)  are  separated 

by  -  uj^  >  0,  and  the  frequency  mismatch  with  the  detuning,  will  be 

indicated  by  i  *  id.  -  (w  -  id  ) .  The  strength  of  the  dipole  coupling  is 
L  e  8 


expressed  in  terms  of  the  Rabi  frequency 
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n  *  dS  EQ  <e|y.e|g>  , 


(2.9) 


which  contains  the  dipole-moment  operator  p  of  the  vibrational  bond.  The 
interaction  Hamiltonian  is  then  -y.E(t). 

The  full  equation  of  motion  can  now  be  summarized  as 


i-frj^  -  [H,p]  ,  p+  «  p  ,  Trp  *  1  , 


(2.10) 


where  p(t)  is  the  density  operator  of  the  atomic  bond  and  the  phonon  field. 
The  Hamiltonian  is  given  explicitly  by 


» ■  »b ♦  v  - 1  vv  +  4> s 


ks 


-  iio.  t  iw^  t 

-  rfin{  |e><g|e  +  e  j  gXe  |  }  , 


(2.11) 


with  the  free  evolution  of  the  adbond  H^  and  the  phonon  field  H  ^  given  by 
Eqs.  (2.6)  and  (2.4),  respectively.  The  coupling  parameter  with  the  phonon¬ 
mode  ks  is  given  by 


f  *i  v1  )1/2 

ks  ^2Mvm(k)J  ?ks‘  --L 


(2.12) 


III.  DRESSED  STATES 

The  interaction  between  the  bare  states  and  the  laser  mode  can  be 

diagonalized  directly.  To  this  end  we  introduce  the  transformed  density 
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operator  oft)  by 


-  iui.  tP  iw.  tP 

o(t)  *  e  8  p(t)  e  8  ,  (3.1) 

with  P^  the  projector  onto  the  ground  state.  This  transformation  eliminates 
the  oscillatory  factors  in  the  interaction  term  -p.E(t),  but  at  the  expense 
of  the  appearance  of  a  time-dependent  coupling  with  the  phonon  field.  The 
equation  of  motion  for  o(t)  reads 


.*  *  f't. 


“g  ■  IHd  +  HPh  •  01 


l  6ks  ((aks  +  aks>  <1<t)  •  ■ 


(3.2) 


-iuj.  tP  iuL  tP 
q(t)  *  e  ^  S  e  ® 


The  dressed-atom  Hamiltonian  H,  is  defined  as 

d 


(3.3) 


\  /fiuJiPi 


i*e,g 


+  ^(tuL  +  ^e  +  t°g)(Pe  +  V  ' 
-  |eXg|  +  |g><e| )  , 


(3.4) 


which  has  the  interpretation  of  the  bare  atom  dressed  with  the  photons  of 

25 

the  laser  field,  including  the  interaction. 


Similarly  to  Eq.  (2.6),  we  now  write 


H .  *  )  'fuL.P.  , 

d  L  11 


(3.5) 


with  fio^  the  eigenvalues  of  H^,  and  the  projectors  onto  the  eigenstates 
|i>  of  H^.  The  eigenvalue  equations  of  are  readily  derived  from  Eq. 
(3.4).  We  obtain 


H.|i>  *  dW|i>  ,  i  =  bare  , 

a  l r 


Hd|t>  =  , 


i  *  e,g 


(3.6) 


(3.7) 


ai.  ■  4(  u),  +  ui  +  (jj)+4n'  , 


(3.8) 


n'  -  a(i  +  cV)1/2 


(3.9) 


Hence  the  dressed  states  are  identical  to  the  bare  states,  provided  that  i  * 
e,g.  The  states  |e>  and  |g>  are  coupled  to  form  |+>  and  |->  according  to 


|+>  ■  |e>  cos0/2  +  U>  sin0/2  , 
l*>  *  U>  cos0/2  -  je>  sin0/2  , 


(3.10) 

(3.11) 


which  is  a  parametrization  in  terms  of  the  angle  0/2,  where  0  is  defined 


0  *  arctan(fl/A) 


(3.12) 


For  weak  driving  fields  (n  -*•  0)  we  have  |+>  -  |e>  and  |->  -  |g>,  but  for 
strong  fields  we  find  |±>  r  (|g>  ±  |e>)//Z.  The  summation  in  Eq.  (3.5)  runs 
over  the  bare  states  i  *  e,g,  and  over  i  *  ±,  where  P+  *  |+><±|.  This 
diagonalizes  the  laser  interaction.  The  relation  between  the  eigenvalues  of 
and  is  illustrated  in  Fig.  1. 

The  coupling  to  the  phonon  field  is  now  regarded  as  an  interaction 

between  the  dressed  bond  and  the  free  phonons.  The  time  dependence  of  this 

2 

interaction  is  embodied  in  q(t)  from  Eq.  (3.3).  With  the  property  P  =  P  , 

g  S 

we  can  expand  the  exponentials  as 


iiui.  tP  tied,  t 

.^K=l-P+e  LP  , 
g  g 


(3.13) 


which  allows  us  to  write 


-  iuL  t  iu).  t 

q(t)  *  S  +  e  LPS+eLSP 


(3.1A) 


Here  we  used  that  <g|S|g>  =  0,  and  we  introduced  the  time-independent  part 
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of  the  interaction  by 


5  *  S  -  P  S  -  SP 
g  g 


(3.15) 


The  identification  of  the  contributions  to  the  interaction  of  the  three 
frequencies  in  Eq.  (3.14)  will  facilitate  the  following  computations 
considerably. 


IV.  INTERACTION  PICTURE 

In  contrast  to  the  photon  field,  where  only  the  laser  mode  is  occupied, 
the  phonon  field  consists  of  many  modes  ks,  which  all  contribute  to  the 
interaction.  This  prohibits  an  immediate  diagonalization,  so  other  methods 
have  to  be  applied.  In  this  section  we  set  up  the  notation  and  write  down 
the  basic  equation.  To  this  end,  we  introduce  the  density  operator  in  the 
interaction  picture  by 


r(H,  +  H  )t  4(H  +  H  )t 

S(t)  -  e*  d  ph  o(t)  e*  d  Ph 


(4.1) 


When  there  is  no  coupling  between  the  dressed  system  and  the  phonon  field, 
we  have  o(t)  =  o(0)  =  o(0).  The  equation  of  motion  for  o(t)  follows  from 
Eq.  (3.2),  and  we  obtain 


ifi  ^  5(t)  *  [W(t)  ,  o(t)) 


The  interaction  Hamiltonian  W(t)  is  then 


(4.2) 


uf .  \  _  Vo  (  *iw(k)t  ,  t  iu)(k)tx  . 

w(t)  "  L  Ks  (aks  e  +  aks  e  }  q(t)  ’ 


with  the  transformed  atomic  part  given  by 


q(  t )  =  e  q( t )  e 


(4.3) 


(4.4) 


If  we  expand  the  exponentials  in  the  projectors  on  the  dressed  states  and 
use  the  expression  (3.14)  for  q(t),  we  obtain 


q(t)  -  l 


id.  t  -iuLt  iwTt 

e  Kt  Pk  (5  +  e  L  PgS  +  e  L  SPg}  ^  . 


(4.5) 


where 


*kl  "  "k  '  “s 


(4.6) 


is  the  level  separation  between  the  dressed  states  |k>  and  |t>.  Besides  the 
usual  time -dependence  exp^&^t)  in  the  interaction  picture,  we  now  also 
find  exp(i(Akl  ±  u^Jt). 

Since  the  system  of  interest  is  the  adbond,  we  define  the  reduced 
density  operator  of  the  dressed  atom  by 


°0(t)  "  TrPh  o(t)  ’ 


(4.7) 


and  similarly  for  Og(t).  Taking  the  trace  over  the  phonon  field  in  Eq. 
(4.1)  yields 


/ '•ftHdt  -  ,  .  h1!dt 
oQ(t)  =  e  c>q( t )  e 


(4.8) 


and  therefore  it  is  sufficient  to  derive  an  equation  for  Og(t).  If  we 
iterate  the  first  integral  of  Eq.  (4.2)  twice,  differentiate  the  result  with 
respect  to  time,  and  take  the  trace  over  the  phonon  field,  we  obtain 


dt  °0^  =  "  Tr-*’  ^ 


2  ph 


dt'  (W(t)  ,  ;«(t')  ,  5(t')]J  . 


(4.9) 


The  right-hand  side  of  this  equation  accounts  for  the  phonon  interaction  in 
the  evolution  of  o„(t). 


V.  PHONON  RESERVOIR 
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The  phonon  field  will  be  considered  as  a  large  reservoir,  or  heat 

bath,  with  a  short  recovery  time.  Then  the  integral  in  Eq.  (4.9)  can  be 
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evaluated  by  standard  techniques,  so  here  we  merely  state  the  result. 

Care  should  be  exercised,  however,  since  the  atomic  part  of  the  interaction, 
q(t),  contains  three  frequencies  for  a  single  kfc  matrix  element.  This  gives 
rise  to  a  number  of  cross  terms,  which  do  not  appear  without  a  laser  field. 

The  average  number  of  phonons  in  mode  w  for  the  crystal  in  thermal 
equilibrium  at  temperature  T  is 

n(u)  =  [exp(fWk  T)  -  1]  *  .  (5.1) 

D 

We  extend  the  definition  of  n(w)  to  negative  values  of  w,  which  enables  us 
to  write 

n(w)  +1  *  -n(-w)  .  (5.2) 


This  will  simplify  the  notation  considerably.  Next  we  introduce  another 
function  of  w  by 

fU)  «  Q  ^  l^djjf^l  k2(M}  [n(u))  *  1]  •  (5-3) 

with 


C  2irM 


(5.4) 


in  terms  of  the  dispersion  relation  k  =  k(a>),  which  is  the  inverse  of  u>  = 
ou(k).  This  f(«)  enters  because  of  the  summation  over  all  modes  ks  of  the 
reservoir . 

After  these  preliminary  definitions,  we  can  evaluate  the  integral  in 
Eq.  (4.9).  We  obtain  the  simple  result 
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dt  °0^  2  \  l5k«J  f^kH^  +  °0^t^k  "2^g,<k  I  °o^  I k>* 

kH 

~  2  \  8±  lSg*l  f^±fc  ‘  “l* 

±11 

*  (f±  50(t)  +  °0(t)?±  -  2fi<±|o()(t)|±>} 

'  1  l  g±  lsgk|2  f(*k±  +  V 

k± 


x  (&k50(t)  +  oQ(t)?k  -  2£+<k|oQ(t) |k>}  . 

-g+g_  lSeg|2(f  (uJL)  +  f('t\))  l  P±°0(t)  P -  •  (5>5) 

±  + 

The  first  term  comes  from  the  time-independent  part  5  of  the  interaction, 

2  2 

and  is  proportional  to  1 5k ^ |  =  |<k|5|H>|  .  Here  the  summation  runs  over 
all  dressed  states  k,«..  The  second  and  the  third  terms  only  connect  the 
dressed  states  ±  with  all  other  states.  In  the  first  term,  the  overall 
factor  contains  f(Ak^),  which  depends  only  on  the  level  distance 
but  in  the  other  two  terms,  this  frequency  is  shifted  by  the  optical 
frequency  Furthermore,  they  involve  the  geometrical  factors 

gt  =  <±|Pgl±>  .  (5.6) 


e.g.,  the  matrix  elements  of  the  projector  on  the  ground  state  with  respect 
to  the  dressed  states.  This  factor  is  only  non-vanishing  for  the  states  |+> 
and  |->,  which  explains  the  limited  summations  in  the  last  two  terms.  For 
later  purposes  we  note  that  g+  can  be  expressed  in  terms  of  the  optical 
parameters  O'  and  A  according  to 


_  A1  +  A 
8±  ~  20' 


(5.7) 


as  can  be  found  from  the  definitions  (3.10)  and  (3.11)  of  the  dressed 


states . 


Notice  that  the  operator  equation  (5.5)  for  Og(t)  involves  only  well- 

defined  quantities.  The  projectors  inside  the  curly  brackets  select  a 

specific  matrix  element  of  Og(t)  with  respect  to  the  dressed  states,  and  the 

prefactors  contain  f(w),  a  given  function  of  ui  and  the  temperature,  and  the 

matrix  elements  of  S  =  dV^/dz.  Besides  that,  the  optical  parameters 

2 

appear,  which  are  known  functions  of  £1  and  the  detuning  A.  In  view  of  Eq. 

2 

(2.9),  the  parameter  ft  is  proportional  to  the  laser  intensity  and  to  the 
square  of  the  transition-dipole  moment.  It  is  this  combination  that 
determines  whether  a  laser  field  can  be  treated  as  a  weak  field 
(perturbation),  or  not.  We  address  this  item  in  more  detail  in  Section  VII. 
The  temperature  enters  parametrically  through  f(u>)  from  Eq.  (5.3).  If  the 
laser  heating  of  the  crystal  is  considerable,  the  temperature  will  depend  on 
time,  and  hence  an  additional  macroscopic  heat-transport  equation  has  to  be 
supplied.  However,  if  the  direct  heating  is  dominant,  the  details  of  the 
atomic  evolution  and  the  radiative  excitation  of  the  adbond  have  no 
importance.  Therefore,  we  assume  an  "adiabatic  following"  with  temperature 
in  this  paper. 


VI.  THE  BLOCH  EQUATIONS 

For  a  given  transition  |k>  -*•  1 1>  there  might  be  a  contribution  from 
more  than  one  term  in  Eq.  (5.5)  if  k  =  t  or  t  =  i.  In  order  to  distinguish 
between  the  various  contributions  in  every  transition,  we  rearrange  the 
terms.  Furthermore,  we  will  express  all  matrix  elements  of  S  and  5  in  terms 
of  bare-state  matrix  elements.  This  implies  that  the  optical  parameters 
will  appear  in  the  overall  factors,  and  hence  the  effect  of  the  laser  is 
tracked  down  explicitly  in  every  transition. 


From  the  definition  of  5,  Eq.  (3.15),  it  follows  immediately  that 


6. 
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so  that  for  states  |k>,|£.>,  which  are  not  directly  coupled  by  the  laser,  we 
can  replace  5  by  S.  If  one  of  these  states  equals  |+>  or  |->,  we  obtain 

|S±kl2  “  8 _ I sek 1 2  *  «  k  *  ±  »  (6.2) 

and  again  g_  accounts  for  the  optical  effects.  We  note  that  these  matrix 

+ 

elements  connect  the  excited  state  |e>  with  |k>.  Finally,  we  find 


<k|5| fc>  =  0 


(6.3) 


if  both  |k>  and  | t>  are  one  of  the  states  |+>  or  |->.  For  the  matrix 
elements  of  S  itself,  as  they  appear  in  the  second  and  the  third  terms  of 
Eq.  (5.5),  we  find 


lsgll2  -  8.|seg|2 


(6.4) 


A  transition  |k>  -»  | t>  might  gain  contributions  from  more  than  one  term 
in  Eq.  (5.5),  but  for  the  time  evolution  of  5g(t)  these  terms  add  up. 

Hence,  we  can  define  the  effective  rate  constants  a^  for  every  specific 
transition  |k>  >  |t>.  The  radiationless  coupling  between  two  bare  states  is 
governed  by 


V  -  f<  V  lsktl2  •  k  * 1 


l  *  ± 


(6.5) 


For  k  *  l  this  rate  constant  is  zero,  because  we  assumed  =  0.  If  one  of 
the  states  equals  |+>  or  |->,  we  obtain 


■it  ■  *-  I8./  *  «,  f(iit  '  “L>  lSgt!2  • 

Si  "  8;  lSeJ2  +  *i  f<Ski  +  V  'V2  ’ 


) 

I 


(6.7) 
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for  t  *  t,  k  *  t.  The  transition  rates  between  the  two  |i>  doublets  are 
determined  by 


*++ 

•  a_  ■ 

8+  8 

_[f(-“»L)  + 

«“l»  Is.,!2 

1 

(6.8) 

V 

-L+2 

f  (-Q* 

2  f(-Q*  +  u^)) 

Is.,'2  • 

(6.9) 

-  Ig_2 

f  (O’ 

‘  V  +  *+ 

f(Q'  +  u^)]  \l 

5e.'2  • 

(6.10) 

We  note 

that  Sk 

*  0 

if  k  -  +  or 

k  - 

With  these 

identifications 

of  the  various 

contributions,  Eq. 

(5.5)  can 

be  summarized  as 

d 

dt 

50(t)  - 

•H 

akt  lPk°0(t)  +  °0(t)Pk  • 

2Pt<k|50(t)|k>] 

kt 

a++ 

l  Pt  5o«> 

P_ 

+ 

(6.11) 

for  the  atomic  density  operator  in  the  interaction  picture.  Transforming 
back  with  Eq.  (4.8)  then  yields  the  optical  Bloch  equations  for  this  system. 
We  find 

itf  ^  o0(t>  ■  [Hd,oQ(t)]  -  iflrOgU)  ,  oQ(t)*  *  oQ(t)  ,  TroQ(t)  -  1  , 

(6.12) 

where  the  Liouville  operator  r  is  defined  by  its  action  on  an  arbitrary 
operator  p  according  to 

rp  "  2  ^  ®kt^kp  +  p^k  *  +  *+_►  ^  p  (6.13) 

kt  i  + 

Here  the  summation  runs  over  all  dressed  states.  This  T  represents  the 
relaxation  of  the  adbond  due  to  single-phonon  transitions,  and  it 
incorporates  the  combined  phonon-photon  processes,  as  will  be  explained  in 


the  next  section. 
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VII.  DIAGRAMS 

The  structure  of  the  relaxation  operator  T  is  quite  transparent.  Every 

kt-combination  assumes  the  form  2^ic°0  +  °0^k  "  2f£<k|og|k>)  t^mes  a 

constant  a,  „ ,  and  in  the  next  section  it  will  become  clear  that  this  action 
kt 

on  Og  corresponds  to  a  transition  jk>  -*■  ji>  at  a  rate  of  times  the 
population  of  level  |k>.  Therefore,  a^  is  termed  the  rate  constant  or  the 
inverse  lifetime  of  the  transition.  On  the  other  hand,  we  actually  know 
which  term  corresponds  to  which  transition,  because  every  term  reflects  a 
single-phonon  excitation  or  decay.  The  particular  transition  can  then  be 
inferred  from  the  definition  (6.5)-(6.10)  of  the  rate  constants,  since  any 
appearance  of  a  factor  f(<n)  displays  the  occurrence  of  a  single-phonon 
transition  with  phonon  frequency  jwj.  Here  a  positive  argument  u>  pertains 
to  a  downward  transition,  and  w  <  0  corresponds  to  an  excitation.  Both  f(u>) 
and  f(-u>)  appear  for  every  transition,  and  the  inverse  lifetime  for  a  decay 
is  always  larger  than  for  an  excitation.  For  zero  temperature  we  have  n(u>) 

*  0,  so  that  the  upward  transitions  disappear  identically,  and  we  are  left 
with  a  pure  decaying  system.  The  occurring  transitions  between  dressed 
states  with  their  optical  part  of  the  rate  constants  are  pictorially 
represented  in  Fig.  2,  and  in  Fig.  3  we  plot  the  explicit  dependence  on  the 
optical  parameters  ft  and  6. 

At  this  stage  it  is  elucidating  to  consider  the  limit  of  a  weak  driving 

2  2 

field.  The  laser  intensity  is  proportional  to  ft  .  For  ft  =  0  we  find  from 
Eq.  (5.7) 


8 


(7.1) 


for  the  parameters  which  determine  the  transition  rates.  Furthermore,  the 
dressed  states  of  Fig.  1  approach  the  bare  states,  since  ft'  -  A.  If  we 


inspect  the  transitions  of  Fig.  2,  we  can  now  also  consider  them  as  bare- 
state  transitions,  where  every  appearance  of  a  factor  g+  corresponds  to  one 
photon.  This  is  illustrated  in  Fig.  4  for  single-phonon  excitations.  We 
obtain  three  diagrams  without  photons,  and  five  radiation-assisted 
transitions,  which  disappear  for  a  zero  intensity. 

VIII.  MATRIX  ELEMENTS 

Equation  (6.12)  contains  all  information  about  the  state  of  the  bounded 
atom,  including  its  interaction  with  the  phonon  field  and  the  radiation.  It 
is  a  condensed  operator  equation,  but  in  order  to  solve  it,  we  have  to 
rewrite  it  in  matrix  form.  This  is  easily  accomplished  if  we  use  the 
explicit  form  (6.13)  of  the  damping  operator.  The  equation  of  motion  for 
the  kt-th  matrix  element  of  oQ(t)  attains  the  form 

^  <k|oQ(t)|l>  -  -  -  it)  <k|oQ(t)|t> 

-  i  <k|,0<t)|i>  £  ♦  .ln> 

n 

+  5kt  \  *nk  <nl°o(t)ln>  • 
n 

-  <k|o0(t>|l>  l  «t-  •  (8'° 

± 

Now  we  can  distinguish  two  cases.  First,  for  k  *  t  and  k  and  t  not  both 
equal  to  +  or  -  we  find 

~  <k | °q( t ) { t>  »  {-i(«k  -  <*>t)  *  |(A^  +  Afc)}  <k|og(t)|l>  ,  (8.2) 

where  we  introduced  the  abbreviation 


n 


Equation  (8.2)  is  an  equation  for  a  single  coherence  of  oQ( t),  and  we  note 
that  its  time  evolution  decouples  from  the  evolution  of  the  other  matrix 


elements.  The  solution  is 


-i(d).  -  w  )t  -t(A.  +  A  )t 

<k|o0(t)|t>  *  e  K  e  Z  K  1  <k|O()(0)|t>  (8.4) 

for  a  given  initial  state  Oq(0).  For  t  >>  (A^  +  A^)  *  we  find  <k. { 0q( t )  |  «.> 
0  for  any  initial  state  Oq(0).  This  justifies  the  random-phase 
approximation  with  respect  to  dressed  states  and  in  the  long-time  limit. 
Since  the  coherences  between  dressed  states  vanish,  they  do  not  disappear 
with  respect  to  the  bare  states,  as  is  commonly  assumed  in  a  random-phase 
approx imat ion . 

For  k  *  l  we  obtain  the  master  equation 

^  "k(t)  *  •  Vk(t)  + 1  atknt(t)  •  (8-5> 

l 

where  we  have  denoted  the  population  of  the  k-th  level  by 

n^(t)  -  <k|oQ(t)|k>  .  (8.6) 

Equation  (8.5)  involves  the  populations  of  all  levels,  but  it  does  not 
couple  with  the  coherences.  This  shows  that  a  master-equation  approach  is 
only  correct  with  respect  to  the  dressed  states,  since  a  transformation  of 
Eq.  (8.5)  to  bare  states  involves  coherences.  In  the  bare-state  basis  we 
have  to  include  all  matrix  elements  simultaneously,  which  is  a  cumbersome 
procedure. 

Equation  (8.5)  for  the  populations  is  a  simple  set  of  linear  first- 
order  differential  equations  with  time- independent  coefficients,  and  it  can 
be  solved  easily  for  any  configuration  of  states  |k>.  The  set  (8.5)  has  to 
be  accompanied  by  the  normalization 


y  nu(t) 


Tro_(t)  -  1 


(8.7) 


The  structure  of  Eq.  (8.5)  is  obvious.  The  population  of  |k>  decays 
exponentially  with  an  inverse  lifetime  A^  due  to  decay/excitation  to  all 
other  levels,  and  it  is  modified  by  the  gain-term  a^^n^(t)  from  every  other 
level  |l>.  This  shows  that  a  |k>  -*■  j t>  transition  occurs  at  a  rate 
*k£nj[(t)»  as  already  mentioned  in  Section  VII.  Furthermore,  level  |k> 
decays  with  a  total  rate  A^n^(t),  and  the  coherence  between  |k>  and  |t> 
decays  with  ^(A^  +  A^).  This  also  clarifies,  at  that  stage,  the  arbitrary 
factor  which  was  separated  from  the  overall  factors  in  Eq.  (5.5). 

IX.  DECAY  OF  DRESSED  STATES 

The  rate  constants  a^g  for  every  transition  |k>  ■*  1 1>  between  dressed 
states  are  defined  in  Eqs.  (6.5)-(6. 10).  They  are  also  defined  for  k  or  l 
equal  to  e  or  g  in  Eq.  (6.5).  We  now  make  a  slight  approximation  in  the 
level-distance  dependence  of  the  constants  by  assuming  that  f(u>)  does  not 
vary  significantly  over  one  |±>-doublet.  This  is  exact  for  weak  fields  and 
perfect  resonance.  Then  we  can  express  the  rate  constants  entirely  in  the 
bare-atom  transition  constants  (6.5)  and  the  laser  parameters  g±.  We  obtain 

*tt  ■  B;  V  +  B,  V  •  ‘  *  *  '  (9'1) 


St  ■  *-  Se  +  Bt  akg 


k  *  t 


(9.2) 


a.  .  ■  a 


g,  g  U  o  +  a  )  , 

+  -  eg  ge 


+- 


8+  age  +  8.  a 


eg 


2  ,  2 
a-+  *  «-  a*e  +  8+  at 


(9.3) 

(9.4) 

(9.5) 


ge  eg 

In  analogy  with  the  definition  (8.3)  for  the  total  decay  constant  A^  of  a 


dressed  state,  we  now  define 


which  would  be  the  relaxation  constant  for  the  bare  state  |k>  if  the  laser 
field  were  not  present.  Now  it  can  be  proven  immediately  from  Eqs.  (9.1)- 
(9.5)  that  the  definition  (9.6)  is  identical  to 


A^  ■  £  a^n  ,  k  *  bare  ,  (9.7) 

n  *  dressed 

provided  that  we  define  a^±  and  a^  as  in  Eq.  (9.2).  This  remarkable  result 
states  that  the  decay  constant  for  the  bare  state  |k>,  including  k  =  e  and  k 
*  g,  is  unaffected  by  the  presence  of  the  laser. 

The  decay  of  the  dressed  states  |±>  can  now  easily  be  expressed  in 
terms  of  the  rate  constants  of  the  bare  levels  |e>  and  |g>.  The  result  is 


8;  *.  +  *.  *8 


(9.8) 


Here  A£  and  A^  are  defined  by  Eq.  (9.6),  and  they  are  independent  of  optical 
parameters.  The  g^-f actors  enter  into  Eq.  (9.8),  so  that  the  decay  of  the 
dressed  states  does  depend  on  0  and  A.  From  g+  +  g_  a  1  it  follows  that 


A  +  A  *  A  +  A  ,  (9.9) 

+  -  e  g 

which  is  not  so  obvious,  since  A  +  A_  or  A  +  A  have  no  direct  physical 
interpretation.  Recall  that  the  loss  rate  from  the  doublet  | ±>  equals 
A+n+(t)  +  A_n_(t). 


X.  STEADY  STATE 

The  master  equation  (8.5)  can  be  solved  immediately  for  any  initial  set 
of  populations  n^(0).  The  transient  behavior  of  the  solution  is  however 
trivial,  since  we  have  a  pure  exponentially-decaying  system.  Furthermore, 


it  is  experimentally  not  feasible  to  prepare  the  system  in  a  well-defined 


initial  state.  Due  to  the  relaxation  however,  every  °q(0)  will  evolve  to  a 
steady-state  o,  defined  as 


o 


lim 

t  4  CD 


°0(t) 


(10.1) 


which  is  not  necessarily  unique.  Similarly,  we  denote  n^®)  by  n^.  In  this 
limit,  the  master  equation  reduces  to 


1 


l  “tk5* 

*  dressed 


(10.2) 


or  equivalently 

l  (st  aik  '  “k  aki>  ’  0  • 

t  *  dressed 


(10.3) 


We  will  now  derive  some  salient  features  of  the  steady-state  solution. 
We  recall  that  the  coherences  vanish  with  respect  to  the  dressed  states,  so 
that  o  is  given  by 


k  *  dressed 

in  terms  of  the  solution  of  Eq.  (10.3)  and  the  projectors  on  the  dressed 
states.  With  the  definition  of  the  |i>  states  in  Eqs.  (3.10)  and  (3.11), 
we  can  transform  the  density  matrix  a  to  the  bare-state  representation, 
which  yields 

o  -  ^  n^  Pk  +  (n+  -  n_)(|e><g|  +  |g><e|)  .  (10.5) 

k  *  bare 

Here  n^  is  the  same  as  n^  in  Eq.  (10.4)  for  k  *  e,g,  and  the  relation 
between  the  populations  n  ,  n  and  n  is  explicitly 
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ne  "  8.n+  +  g+n_  , 


(10.6) 


ng  "  8+n+  +  g_n_ 


(10.7) 


Equation  (10.5)  shows  that  a  non-vanishing  coherence  <e|o|g>  appears  in  the 
bare-state  representation.  With  the  aid  of  Eqs.  (10.6)  and  (10.7),  we 
find  the  relation 


n+  -  n 


n  -  n 

_8 _ S 


8+  '  8. 


(10.8) 


and  with  the  expression  (5.7)  for  g±,  we  obtain  for  the  coherence 


<eHg>  =  U  (ne  ’  ng>  * 


(10.9) 


In  the  absence  of  a  laser  we  have  0*0,  and  <e  1 0 1 8>  vanishes . 


Relations  (10.6)  and  (10.7)  can  be  inverted  to  express  n^  in  terms  of 


n@  and  n^.  If  we  subsequently  use  Eq.  (9.8)  for  the  relation  between  the 


rate  constants,  we  find 


A, n,  +An  *  A  n  +An 
+  +  --  ee  gg 


(10.10) 


This  identity  expresses  that  the  total  loss  rate  of  |±>  equals  the  sum  of 
the  loss  rates  of  |e>  and  |g>,  as  could  be  expected. 

An  advantage  of  the  application  of  dressed  states  is  that  the  master 
equation  (10.3)  has  a  simple  form.  With  the  relations  of  Section  IX,  we  can 
transform  Eq.  (10.3)  to  its  bare-state  equivalent.  We  find 
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where  all  quantities  pertain  to  bare-state  properties.  This  result  is 

remarkably  simple  also.  The  laser  only  modifies  the  master  equation  for  k  * 

g  and  k  *  e  in  the  form  of  an  additional  term,  and  the  optical  parameters 

2  2 

only  appear  in  the  combination  Q  /A  .  The  set  (10.11)  couples,  however, 
between  all  k-values,  and  therefore  the  extra  terms  affect  all  populations 
iv^.  We  notice  that  the  right-hand  side  of  the  master  equation  for  k  *  e  and 
k  *  g  contains  the  same  radiative  contribution.  If  we  add  the  two 
equations,  we  are  left  with 


)  n.  (a.  +a.)=An+An 
L  *  *g  Ke  g  g  e  e 

k  *  bare 


(10.12) 


and  the  laser  has  disappeared  completely  in  this  combination. 

Let  us  now  consider  a  specific  situation.  Suppose  that  the  ground 
state  |g>  is  the  lowest  state  and  that  its  separation  from  every  other  level 
is  larger  than  the  Debye  frequency.  Then  the  single-phonon  transition  rates 
from  and  to  |g>  can  be  neglected  in  comparison  with  the  radiative 
coupling  to  |e>.  From  Eq.  (10.11)  for  k  ■  g,  it  then  follows  that 


n  -  n  ,  (10.13) 

g  e 

and  the  master  equation  reduces  to 

l  Vtk  ■  Vk  •  <10-14) 

l  -  bare 

for  every  bare  state  |k>.  This  is,  however,  the  Master  equation  in  the 

absence  of  the  laser,  and  it  might  seem  that  the  laser  has  no  effect  at  all. 

Such  is  not  the  case.  The  master  equation  always  provides  one  superfluous 

equation,  which  allows  us  to  impose  the  normalization  (8.7)  on  the  solution. 

Now  we  have  the  restriction  n  =  n  ,  so  Eq.  (8.7)  should  be  replaced  by  the 

e  g 
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This  proves  that  for  this  specific  configuration  the  distribution  of 
population  over  the  states  is  purely  thermal  for  k  *  g,  e.g.,  independent  of 
the  laser  intensity.  The  normalization  (10.15)  indicates  that  we  can  ignore 
the  ground  state  completely,  provided  that  we  take  into  account  that  there 
is  a  population  equal  to  ng  missing.  The  desorption  rate  is  determined  by 
the  populations  of  the  upper  levels,  which  implies  that  the  presence  of  the 
laser  enhances  the  desorption.  Indeed,  population  which  would  be  trapped  in 
the  ground  state  is  now  continuously  pumped  upwards  to  the  desorptive 
states . 


XI.  SIMPLE  THREE -LEVEL  EXAMPLE 

In  order  to  demonstrate  the  applicability  of  our  approach,  we  elaborate 
an  example.  Suppose  that  we  have,  beside  the  levels  |e>  and  |g>,  one  other 
level  |n>.  Three  independent  rate  equations  are  then 
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which  are  Eq.  (10.11)  for  k  =  g,  k  *  n  and  the  normalization.  If  we  add  Eq 
(11.1)  to  Eq.  (10.11)  for  k  =  e,  we  obtain  Eq.  (11.2).  The  solution  of  the 

set  ( 1 1 . 1 ) - ( 1 1 . 3 )  is  easily  found.  For  n  we  obtain 
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and  similar  expressions  hold  for  n^  and  n^. 
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the  dependence  on  the  optical  parameter  fi  /& 
inverse  lifetimes. 

If  we  now  set  u  -  u  >  ul  and  ui  -  u  > 

n  g  O  eg 


This  result  displays  explicitly 
and  the  laser- independent 

ujp,  Eq.  (11.4)  reduces  to 
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and  the  populations  of  ng  and  n^  are 
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This  shows  that  in  this  limit  the  dependence  on  the  laser  intensity  indeed 
vanishes,  as  stated  in  general  in  the  previous  section.  Without  a  laser, 
the  stationary  state  for  this  limit  becomes 


n  ■  n  *  0  ,  n  «  1  .  (11.7) 

n  e  g 


Hence  the  laser  sustains  a  finite  population  of  the  states  |e>  and  |n>, 
which  have  a  finite  desorption  rate. 


XII.  CONCLUSIONS 

We  have  studied  the  interaction  of  a  physisorbed  atom  with  strong 
coherent  radiation.  The  mechanisms  of  spontaneous  phonon  transitions  of  the 
vibrational  atomic  bond  and  the  photon  absorptions  and  stimulated  emissions 
all  give  rise  to  a  redistribution  of  the  level  populations.  The  thermal 
relaxation,  due  to  phonon-exchange  reactions,  drives  the  system  towards 
equilibrium  with  the  phonon  reservoir,  but  the  irradiance  tends  to  maintain 
a  temperature  gradient  between  the  adbond  and  the  phonon  field.  This  can  be 
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I 

understood  from  the  underlying  photon-phonon  exchange  reactions  via  the  ! 

i 

adbond,  as  pictorially  displayed  in  Fig.  4.  Since  the  absorption  rate  does  j 

not  necessarily  equal  the  rate  of  stimulated  emission,  a  net  photon 
absorption  rate  will  remain.  For  strong  fields  this  simple  picture  does  not 
apply  anymore,  and  we  have  to  consider  multiphoton  processes  as  well.  This  | 

was  accomplished  by  the  introduction  of  the  dressed  states,  which  are 

t 

schematically  presented  in  Fig.  1.  The  Hamiltonian  of  the  atom  plus  the 
laser,  including  the  interaction,  is  diagonal  with  respect  to  these  dressed  j 

states.  The  interaction  with  the  phonon  field  was  restricted  to  single¬ 
phonon  transitions,  which  were  assumed  to  be  dominant.  For  situations  where 
this  is  not  the  case,  higher-order  processes  can  easily  be  included.  In  j 

practical  situations  only  two,  not  necessarily  adjacent,  levels  will  be 
coupled  resonantly  by  the  laser.  The  generic  idea  is  to  drive  a  jg>  -*•  |e> 
transition,  with  |g>  as  one  of  the  low-lying  states  and  |e>  as  close  as  j 

possible  to  the  continuum.  This  configuration  prohibits  resonant  coupling 
of  other  transitions  and  gives  a  maximum  enhancement  of  the  desorption. 

Besides  that,  the  optical  frequency  will  be  larger  than  uj^,  and  consequently  I 

the  direct  crystal-heating  is  absent.  Hence,  our  two-level  and  adiabatic 
approach  can  be  regarded  as  quite  general. 

The  dressed  states  appear  to  be  a  convenient  basis  set  for  the  I 

derivation  of  the  optical  Bloch  equations  for  this  configuration.  The 
master  equation  was  obtained  by  taking  the  matrix  elements  of  the  full 
operator  equation,  and  it  was  shown  that  the  populations  obey  the  Pauli 
equation.  The  coefficients  involve,  however,  the  optical  parameters  Q  and 
A,  and  hence  the  lifetimes  of  the  transitions  depend  on  the  properties  of 
the  driving  field.  Furthermore,  it  was  shown  that  the  coherences  with 
respect  to  the  dressed  states  decay  exponentially  to  zero.  The  master 


equation  in  the  bare-state  representation  and  in  the  steady  state  was 
obtained  by  applying  the  relation  between  dressed  states  and  bare  states. 

We  found  that  the  master  equation  greatly  resembles  the  radiationless 
equations,  and  that  the  presence  of  the  laser  could  be  incorporated  entirely 
by  the  addition  of  a  few  terms. 

The  general  theory  was  exemplified  by  the  simple  case  of  a  three- level 
system.  We  achieved  the  explicit  steady- state  solution  for  the  bare-state 
populations,  and  it  was  pointed  out  that  the  system  can  be  desorptive  due  to 
the  continuous  optical  excitation,  whereas  the  bond  would  be  completely  in 
its  lower,  non-desorbing  state  in  the  absence  of  the  irradiance. 
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Fig.  1.  The  energy  levels  of  the  bare  system  plus  the  free  laser  field  are 
represented  by  the  diagram  on  the  left-hand  side.  The  detuning  A  is  taken 
positive.  The  number  of  photons  in  the  laser  mode  is  indicated  by  n.  If 
the  dipole  coupling  -y.E  is  included,  a  new  diagonalization  yields  the 
dressed  energy-level  system  on  the  right-hand  side.  The  level  separations  A 
and  in  the  bare  system  now  become  fl*  and  where  is  always  larger 
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than  A  due  to  Q  *  0.  This  is  the  ac-Stark  shift.  The  state  |k>,  which 


is  not  coupled  to  another  level  by  the  laser,  remains  unaltered.  In 
general,  we  suppress  the  quantum  number  n  and  represent  the  dressed  atom  as 


an  infinite  ladder  of  states  |±>,  which  is  the  famous  Jaynes -Cummings  model 

28 


for  the  interaction  of  a  two-level  system  with  a  strong  radiation  field.' 


Fig.  2.  Illustration  of  the  single-phonon  transitions,  whenever  a  (+>  or 
|->  state  is  involved.  The  parameters  near  the  arrows  are  the  optical  parts 
of  the  prefactors  as  they  occur  in  the  transition  rates.  Without  a  driving 
laser,  the  g+  vanishes,  and  hence  only  the  double-arrowed  transitions 
survive.  For  weak  incident  fields,  the  single  arrows  appear,  and  only  for 
sufficiently  strong  fields  will  the  dashed  transition  |+>  -*•  |->  be  present. 
We  note  that  there  are  no  transitions  in  a  single  doublet.  Furthermore,  the 
corresponding  downward  transitions  have  the  same  optical  factor  but  a 
different  f(u>),  e.g.,  a  different  temperature  dependence.  Besides  these 
transitions,  we  have  single-phonon  processes  between  any  two  bare  states  |k> 
and  |l>,  which  are  obviously  independent  of  the  laser  parameters. 

2  2 

Fig.  3.  The  optical  factors  in  the  transition  rates  as  a  function  of  Q  /A  . 

For  strong  fields  or  small  detunings,  both  g+  and  g  tend  to  their 

2  2 

asymptotic  value  of  1/2,  which  is  the  upper  dotted  line,  and  g,  ,  g  and 
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g+g_  will  approach  1/4  in  this  limit.  For  weak  fields  or  large  detunings, 
2  2 

we  have  g_  *  g_  ■  1  and  g+  *  g+  ■  g+g_  M  0.  For  increasing  field 
2 

strengths,  the  g_  transition  will  diminish  very  strongly,  in  favor  of  the 
other  three  processes  between  the  |±>  doublets. 


Fig.  4.  Single-phonon  excitations  in  the  low-intensity  limit.  For  n 
small,  we  can  consider  the  transitions  between  the  dressed  states,  as  given 
in  Fig.  2,  as  photon-assisted  transitions  between  the  bare  states.  The 
diagrams  in  this  picture  have  the  same  order  as  the  transitions  in  Fig.  2. 
Here  every  diagram  connects  two  of  the  three  bare  states  |k>,  |e>  or  |g>, 
since  we  cannot  populate  a  virtual  level  (dotted  lines).  A  straight  line 
indicates  a  phonon,  and  a  wiggly  arrow  is  a  photon  with  frequency  u^.  The 
three  double-arrow  transitions  in  Fig.  2  are  the  three  radiationless 
transitions,  so  they  persist  even  if  we  switch  off  the  laser  field.  The 
single  arrows  from  Fig.  2.  appear  to  correspond  to  one-photon  transitions  in 
this  diagram,  so  they  arise  already  for  low  intensities.  The  diagram  on  the 
right-hand  side  is  a  two-photon  process,  which  consequently  appears  for 
relatively  strong  fields  only.  It  is  the  dashed  g+  transition  form  |+>  to 
|->  in  Fig.  2.  The  phonon-decay  processes  correspond  to  the  inverse 
diagrams  (reverse  the  direction  of  all  arrows).  In  this  figure  we  only  draw 
the  phonon  excitations,  which  correspond  to  the  absorption  of  a  phonon  by 
the  adbond  from  the  crystal.  The  third  diagram  represents  the  absorption  of 
a  photon  from  the  laser  field,  but  all  other  diagrams  give  rise  to 
stimulated  emissions  of  photons. 
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